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ABSTRACT 

The most luminous Spitzer point sources in the 30 Doradus triggered second gener- 
ation are investigated coherently in the 3-8 jum region. Remarkable diversity and com- 
plexity in their natures are revealed. Some are also among the brightest JHK sources, 
while others are not. Several of them are multiple when examined at higher angular 
resolutions with HST NICMOS and WFPC2/WFC3 as available, or with VISTA/VMC 
otherwise. One is a dusty compact H II region near the far northwestern edge of the 
complex, containing a half dozen bright /-band sources. Three others appear closely 
associated with luminous WN stars and causal connections are suggested. Some are in 
the heads of dust pillars oriented toward R136, as previously discussed from the NIC- 
MOS data. One resides in a compact cluster of much fainter sources, while another 
appears monolithic at the highest resolutions. Surprisingly, one is the brighter of the 
two extended "mystery spots" associated with Knot 2 of Walborn et al. Masses are de- 
rived from YSO models for unresolved sources and lie in the 10-30 Mq range. Further 
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analysis of the IR sources in this unique region will advance understanding of triggered 
massive star formation, perhaps in some unexpected and unprecedented ways. 

Subject headings: stars: massive — stars: pre-main sequence — stars: Wolf-Rayet — 
ISM: individual objects: 30 Doradus — Magellanic Clouds — infrared: stars 

1. Introduction 

The first evidence of discrete, current star formation in 30 Doradus (the largest H II region in 
the Large Magellanic Cloud and indeed in the Local Group) was found by Hyland et al. (1992), who 
reported four luminous protostars, and by Rubio et al. (1992), who discovered 17 near-IR sources 
arrayed in arcs to the west and north of the central cluster R136. Prior to that, 30 Dor had been 
considered an evolved H II region without major ongoing star formation, although "extended" 
(unresolved) far IR emission had been reported west and northeast of R136 already by Werner 
et al. (1978). More extensive subsequent investigations with increasingly powerful instrumentation 
have established that in fact, 30 Doradus is a "two-stage starburst" (Walborn & Parker 1992), 
with substantial triggered star formation primarily at the interfaces between the expanding central 
cavity and remnant giant molecular clouds to the west and northeast (Walborn &; Blades 1997; 
Rubio et al. 1998; Johansson et al. 1998; Scowen et al. 1998; Walborn et al. 1999, 2002; Brandner 
et al. 2001). 

The Spitzer Space Telescope has provided the best, unimpeded imaging capability so far at 
longer IR wavelengths than those readily observed from the ground (or from the Hubble Space 
Telescope). Extensive lists of sources in the Large Magellanic Cloud observed by the SAGE program 
(Meixner et al. 2006, see below), including a number within 30 Doradus, have been provided by 
Whitney et al. (2008) and Gruendl &; Chu (2009). Specific analysis of 30 Dor with the Spitzer data 
has concentrated on the diffuse emission (Indebetouw et al. 2009). Here we undertake the first 
detailed study of the individual and global characteristics of the brightest Spitzer point sources in 
30 Dor, in the context of previous work at shorter wavelengths. 

2. Data 

Our study is based upon the four Spitzer /IRAC bands covering the 3.6-8.0 fim region, in 
which the angular resolution ranges from 1'.'7 through 2'.'0, respectively. Because of the endemic 
multiplicity of massive stars and the distance of the LMC (at which 1'' subtends 50,000 AU), the 
highest available resolutions are essential. Thus, we also discuss ^ O'.'l HST NICM0S/WFPC2 
(Walborn et al. 1999, 2002; Brandner et al. 2001) and WFC3 (E. Sabbi et al. in prep.) data, as well 
as - V' groundbased Rubio et al. (1998), 2MASS (Skrutskie et al. 2006), IRSF (Kato et al. 2007), 
and VISTA/VMC (Cioni et al. 2011) data in comparison with Spitzer. Some further characteristics 
of the surveys and data used in our analysis follow. 
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2.1. Spitzer/^AG^ 

The Spitzer Space Telescope Legacy Program "Surveying the Agents of a Galaxy's Evolution" 
(SAGE; Meixner et al. 2006) imaged a 7° x 7° region in the LMC with ah the IRAC (3.6, 4.5, 5.8, 
and 8.0 /xm; Fazio et al. 2004) and MIPS (24, 70, and 160 /xm; Rieke et al. 2004) bands. The MIPS 
data are not used here because of their lower angular resolution. We use the IRAC "Single Frame 
+ Mosaic Photometry" Archive^ 



2.2. VISTA/VMC 

The "VISTA Magehanic Clouds" survey (VMC; Cioni et al. 2011) is a European Southern Ob- 
servatory (ESO) public survey, which consists of deep imaging in the three NIR bands Y (1.02 /xm), 
J (1.25 /im), and Kg (2.15 /xm), covering both Magellanic Clouds and surrounding areas, with the 
VIRCAM (VISTA IR Camera; Dalton et al. 2006) at the ESO VISTA telescope. The first release 
of VMC data includes images of the 30 Doradus Region (a "tile" of about 1?0 x 1?3). 

We retrieved the 30 Doradus images through the ESO Data Archive Facility. For our analysis, 
we used the "pawprint" images (six of which form a "tile"). The released VMC data have been 
processed through the dedicated pipehne of the VISTA Data Flow System (Irwin et al. 2004, 
Emerson et al. 2004). In addition to the images, the VMC survey contains source catalogues, but 
the aperture photometry they provide includes only a few of our sources and may be unreliable in 
regions of bright, variable nebulosity. 

Accordingly, we performed point-spread-function (PSF) photometry of our sources in the VMC 
images, using DAOPHOT (Stetson 1987) software layered in IRAF|H| Stars were detected at a 
minimum 5a level above the background, with special care to isolate faint companions to the 
targets. 

The VISTA photometric system is tied to the Two Micron Ah Sky Survey (2MASS; Skrutskie 
et al. 2006). To derive the photometric cross-cahbrations, we have selected 440 2MASS catalogue 
stars with the lowest photometric uncertainties (<10% in all three bands) and located in areas 
of faint nebulosity, matching them to our VISTA photometric catalogue. The transformation 
equations obtained are 

-^2MASS - JviSTA = 0.06 X ( JyiSTA - ^VISTa) 

and 

^2MASS - ^VISTA = -0.003 X ( JyiSTA " ^VISTa) • 



^See "The SAGE Data Products Description" at jhttp://data.spitzer.caltech.edu/popular/sage/20090922_enhanced/| 
documents/. 

^IRAF is distributed by the National Optical Astronomy Observatory. 
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The coefficients of these equations are very similar to those found by Cioni et al. (2011). Our 
transformations are vahd for the color range of the 2MASS sources used to calculate them, i.e., 
—0.1 < {J — Kg) < 4.2. Table la shows that the 30 Dor sources with both J and Kg magnitudes are 
in the valid range. For sources with only limits on J, the uncertainties in the Kg transformation 
are also very small, due to the small value of the color-term coefficient. (The Y magnitudes in 
Table la remain in the VISTA photometric system and are not used in our analysis below.) 

2.3. IRSF 

The Infrared Survey Facility (IRSF) Magellanic Clouds Point Source Survey (Kato et al. 2007) 
covered a 40 deg^ area in the LMC in J (1.25 /im), H (1.63 /im), and (2.14 /im) bands, with 
the SIRIUS camera at the South African Astronomical Observatory (SAAO). The average seeing 
throughout the survey was 1'.'2 (FWHM). 

2.4. Results 

Figures 1 and 2 display Spitzer and VISTA images of the field investigated at the same scale, 
with the ten brightest Spitzer sources and some other relevant objects identified. It is seen that 
all of them lie within or near the Tarantula filaments to the west and northeast of R136, except 
for the very luminous and previously little known source SI to the northwest. (The hemispheric 
structure of two-stage starbursts, as discussed by Walborn 2002, is well seen.) The photometric 
data taken from the literature, except for VISTA/VMC as discussed above, for these and other IR 
sources of interest are listed in Table la. The alternate designations are from Hyland et al. (1992) 
or Rubio et al. (1998), while the corresponding SAGE and IRSF identification numbers are given in 
Table lb. Figure 3 shows that the 30 Doradus sources are among the most luminous YSOs in the 
LMC, as listed in the caption references. The 30 Dor sources are discussed individually in detail in 
the following section, emphasizing their remarkable structures and diversity. 




Fig. 1. — A composite of the full 30 Doradus Nebula field in two Spitzer and a VISTA/VMC bands. 
North is up and east to the left. The top ten IR sources are labeled in yellow (S), while bright, 
isolated WN stars are in white (R for Radcliffe, V for VFTS). R136 is the bright blue object below 
center, while the older cluster Hodge 301 with late-type supergiants (green) is labeled. Note the 
color range of the IR sources with this broad baseline. 
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Fig. 2. — Same as Fig. 1 but entirely in VISTA/VMC bands, which have higher resolution. Note 
the multiplicity of several IR sources as discussed in the text, as well as the reversals of relative 
magnitudes in the S3 and SIO systems with respect to Fig. 1. 
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Fig. 3. — A Spitzer CMD comparing the brightest 30 Dor sources (red squares) with the LMC 
YSO candidates of Whitney et al. (2008) and Gruendl & Chu (2009), and the confirmed YSOs of 
Shimonishi et al. 2008, Ohveira et al. 2009, Scale et al. 2009, and Woods et al. 2011. The present 
sources are in the brighter region of the distribution, although several of them are multiple as 
discussed in the text. 



Table la. Photometry 



Spitzer 


Alternate 


R.A.(J2000) 


Decl.(J2000) 






VISTA 






IRSF 








IRAC 






Source 


IDs 


m 


(° ' ") 




Y 


J 


Ks 


J 


H 






[3.6] 


[4.5] 


[5.8] 


[8.0] 


SI 




05:38:31.62 


-69:02:14.6 


16 


68(0.1) 


16.10(0.09) 


14.59(0.08) 








10 


21(0.06) 


9.37(0.06) 


7.72(0.06) 


5.73(0.06) 


S2 


IRS W- 11 


05:38:33.09 


-69:06:11.7 






18.80(0.4) 


14.98(0.02) 




16.95(0.09) 


15.34(0.08) 


12 


(0.10) 


10.56(0.07) 


8.95(0.07) 


7.60(0.07) 


S3 




05:38:34.05 


-69:04:52.2 






>19 


14.32(0.03) 




16.41(0.05) 




10 


5 (0.05) 


9.08(0.04) 


7.83(0.04) 


5.88(0.04) 


S3K ... 




05:38:34.69 


-69:04:50.0 


14 


27(0.01) 


13.37(0.01) 


11.71(0.01) 


13.21(0.02) 


12.17(0.01) 


11.69(0.01) 


11 


27(0.05) 


11.29(0.07) 


10.02(0.21) 




S4 


IRSW-30;P3 


05:38:34.60 


-69:05:56.8 


16 


01(0.07) 


15.23(0.03) 


11.32(0.01) 


14.57(0.07) 


13.18(0.03) 


11.32(0.02) 


8 


78(0.03) 


7.80(0.03) 


6.86(0.02) 


5.77(0.04) 


S5 


IRSW-127 


05:38:39.68 


-69:05:37.9 






>19 


13.39(0.05) 






14.0 (0.03) 


9 


47(0.05) 


8.04(0.04) 


6.95(0.02) 


6.0 (0.09) 


S6 


P2 


05:38:41.36 


-69:03:54.0 


17 


34(0.02) 


15.72(0.01) 


11.52(0.01) 


16.25(0.03) 


13.66(0.01) 


11.5 (0.01) 


9 


13(0.03) 


8.39(0.03) 


7.65(0.02) 


6.89(0.02) 


S7A ... 


IRSN-122 


05:38:46.84 


-69:05:05.4 


19 


50(0.25) 


16.83(0.08) 


13.13(0.02) 


16.13(0.09) 


14.81(0.05) 


12.99(0.02) 


10 


26(0.06) 


9.31(0.05) 


8.25(0.05) 




STB ... 


IRSN-126;P1 


05:38:47.02 


-69:05:01.7 






17.59(0.04) 


11.97(0.01) 


15.59(0.1) 


13.75(0.04) 




8 


81(0.04) 


7.84(0.03) 


7.01(0.03) 


6.0 (0.06) 


S8 




05:38:48.17 


-69:04:11.7 
















10 


11(0.06) 


9.33(0.06) 


7.8 (0.06) 


6.01(0.07) 




IRSN-134 


05:38:48.00 


-69:04:12.7 


17 


46(0.10) 


16.68(0.08) 


14.69(0.10) 


16.1 (0.1) 


15.64(0.12) 


14.57(0.11) 














IRSN-137;P4 


05:38:48.20 


-69:04:11.4 


17 


14(0.09) 


16.36(0.05) 


13.41(0.02) 


16.03(0.07) 


14.86(0.05) 


13.28(0.04) 












S9 


IRSN-152 


05:38:49.27 


-69:04:44.4 






>19 


14.51(0.05) 




16.0 (0.09) 


14.14(0.06) 


10 


3 (0.07) 


9.12(0.06) 


7.55(0.08) 


5.96(0.11) 


SIO .... 




05:38:56.58 


-69:04:17.3 
















11 


59(0.09) 




7.98(0.05) 


6.86(0.10) 




SlOA 


05:38:56.31 


-69:04:16.1 


19 


9 (0.2) 


18.68(0.1) 


15.41(0.04) 




16.9 (0.05) 


14.99(0.03) 














SlOB 


05:38:56.60 


-69:04:17.5 






>19 


15.84(0.02) 




18.64(0.26) 


15.81(0.05) 














SIOC 


05:38:56.84 


-69:04:17.0 






>19 


15.31(0.03) 




17.28(0.06) 


15.23(0.03) 












SIOK .. 




05:38:58.38 


-69:04:21.6 


19 


61(0.1) 


18.58(0.04) 


14.47(0.01) 


17.72(0.07) 


16.53(0.03) 


14.55(0.02) 


11 


82(0.05) 


10.77(0.04) 


9.79(0.16) 




Sll .... 




05:38:52.67 


-69:04:37.5 
















10 


99(0.07) 


10.5 (0.07) 


8.6 (0.07) 


6.98(0.08) 




IRSN-169 


05:38:52.72 


-69:04:37.5 


17 


74(0.05) 


16.70(0.03) 


13.68(0.01) 


16.47(0.04) 


15.17(0.03) 


13.84(0.02) 














IRSN-170 


05:38:53.06 


-69:04:37.4 






18.45(0.08) 


15.41(0.05) 


18.12(0.1) 


16.73(0.07) 


15.28(0.07) 
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Table lb. Source Identifications in the Spitzer SAGE and IRSF Surveys 



Spitzer 


Alternate 


SAGE 


IRSF 


Source 


ID'S 


IRAC^ 




SI .... 




053831.62-690214.6 




S2 .... 


IRS W- 11 


SSTISAGEMA J053833.09-690611.7 


05383309-6906118 


S3 .... 




SSTISAGEMA J053834.05-690452.2 


05383407-6904524 


S3K .. 




SSTISAGEMA J053834.69-690450.0 


05383475-6904504 


S4 .... 


IRSW-30;P3 


SSTISAGEMA J053834.60-690556.8 


05383463-6905569 


S5 .... 


IRSW-127 


SSTISAGEMA J053839.68-690537.9 


05383970-6905382 


S6 .... 


P2 


SSTISAGEMA J053841.36-690354.0 


05384138-6903541 


S7A .. 


IRSN-122 


SSTISAGEMA J053846. 84-690505. 4 


05384685-6905058 


S7B .. 


IRSN-126;P1 


SSTISAGEMA J053847.02-690501.7 


05384701-6905019 


S8 .... 




053848.17-690411.7 


multiple 




IRSN-134 




05384800-6904127 




IRSN-137;P4 




05384820-6904114 


S9 .... 


IRSN-152 


053849.27-690444.4 


05384930-6904444 


SIO ... 




SSTISAGEMA J053856.58-690417.3 


multiple 




SlOA 




05385631-6904161 




SlOB 




05385660-6904175 




SIOC 




05385684-6904170 


SIOK . 




SSTISAGEMA J053858.38-690421.6 


05385837-6904214 


Sll ... 




053852.67-690437.5 


multiple 




IRSN-169 




05385272-6904375 




IRSN-170 




05385306-6904374 



^IRAG designations are from the SAGE IRAG SMP Archive (preceded by "SSTIS- 
AGEMA J") or from Gruendl & Ghu (2009). 
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3. Individual Sources 

Figures 4a, b, c display enlargements of the individual source fields, which show a number of 
significant local structures and relationships that are discussed in this section. 

3.1. SI 

We were surprised to find one of the most luminous Spitzer sources embedded in extensive 
diffuse emission at the northwestern extreme of the field, evidently unrelated to R136. SI appears 
to be associated with a larger CO feature catalogued as 30Dor-06 by Johansson et al. (1998), 
and with an H II region listed as No. 889 by Kastner et al. (2008). In the IR images it appears 
as a somewhat diffuse double source, with the fainter barely resolved component toward the SE. 
Fortunately, a recent HST/WFC3 F775W image (GO 12499, PI D.J. Lennon) reveals what lurks 
within (Figure 5): the object appears as a dusty H II knot, here named "The Skull Nebula," 
containing a compact cluster dominated by four / ~ 19-21 stars in an EW line (the brighter IR 
"component") and two more toward the SE (the fainter IR "component") (E. Sabbi, priv. comm.). 
Two more comparable stars within the outer nebulosity to the N and S may be additional members 
or field objects. 

With reference to Figures 1 and 2, we propose an origin of star formation at this location 
triggered by adjacent SN events. To the NE of SI are the three luminous, "isolated" WN stars 
R144, R146, and R147 (Feast et al. 1960), apparently associated with the large kinematical (Chu 
& Kennicutt 1994) and X-ray (Townsley et al. 2006) shell at the north of 30 Dor. This object 
could have been produced by the explosion of another former WR star. To the SW of SI lies the 
20 Myr old cluster Hodge 301 that has produced 40 SN (Grebel &; Chu 2000) and is surrounded by 
another giant shell. Numerous dark globules can also be seen along the interface between the two 
shells in optical images, suggesting that further star formation is likely to occur in these compressed 
filaments — a novel feature of 30 Doradus. 

3.2. S2 

This relatively fainter source is ISRW-11 of Rubio et al. (1998). (It is also NICOTa of Brandner 
et al. 2001, in which the cross-identification with Rubio et al. was missed.) It is located at the 
northern end of a diffuse structure that is bright in the Spitzer image (Fig. 1) but dark in JHK 
(Fig. 2). As noted by Walborn et al. (1999), with appropriate stretch and scale the object appears 
multiple in the NICMOS data, and analysis of its PSF indicates that it is extended, with a Gaussian 
sigma of 1.29 pix vs. 1.04 pix for a point source in K. Here we analyze the S2 photometry as if 
from a single source, with the caveat that it may actually be a tight multiple system. 
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Fig. 4a. — Larger scale cutouts of the present IR sources in two VISTA/VMC and two Spitzer /IRAC 
bands, which reveal additional structural details. North is up and east to the left in each panel. 
Among other things, note here the reversal of relative magnitudes in the S3 system between the 
two wavelength regimes, and the invisibility of S2, S3, and S5 in J. R135 is a luminous WN star. 
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Fig. 4b. — As Fig. 4a. The object "2!' NW of S6 is an LPV that is brighter than the YSO at shorter 
wavelengths, except when in dust-formation episodes. See B.ST images of S8 in Walborn et al. 
(2002) that show the associated cluster of faint IR sources. 

3.3. S3 

This luminous source apparently within a Tarantula filament is at the apex of an inverted 
triangular asterism with the "isolated" WN star R135 (Feast et al. 1960) and another source 
brighter at shorter wavelengths that we designate S3K. The reversal of their relative magnitudes 
between Figures 1 and 2 is striking; see also Figure 4a. We suggest that the massive wind of the 
WN star may have triggered the formation of S3 in the adjacent nebular filament. As further 
analyzed below, S3K may be a star rather than a YSO. 
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Fig. 4c. — As Fig. 4a. See Walborn et al. (1999) for the remarkable extended structure of S9 in 
HST/NICMOS images; VFTS 621 is a very young 02 star. The positions of two water masers 
reported by Ehingsen et al. (2010) in this field are also marked. VFTS 682 is a luminous WN star. 
Rubio et al. (1998) IRSN-169 is discussed with SIO in Section 3.10. 

3.4. S4 

S4 is one of the most luminous sources in 30 Dor at all NIK wavelengths, also IRSW-30 of Rubio 
et al. (1998), NIC03a of Brandner et al. (2001), and P3 of Hyland et al. (1992). As extensively 
discussed based on HST optical and IR images by Walborn et al. (1999, 2002), it is located in 
the head of a large, bright-rimmed dust pillar oriented toward R136. They measured magnitude 
20.3 in V and 19.3 in /; here we report 11.3 in K and 5.8 at 8 fim (Table la). This is a prime 
object for astrophysical investigation, including spectroscopy from visual through IR. IRSW-26 is 
a companion to the SW, but it is more than 3 mag fainter in K and evidently not significant in the 
Spitzer bands (Figure 4a), so we analyze S4 as a single object. 
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Fig. 5.— SI, "The Skull Nebula," as seen in the i/ST/ WFC3/F775W + groundbased compos- 
ite image presented by E. Sabbi in STScI news release 2012-35 ( |http://hubblesite.org/ newscen- 
ter/archive/releases/2012/35/). See text for discussion. 



3.5. S5 



S5 or Rubio et al. (1998) IRSW-127 is located within a prominent boomerang- (or stapler-) 
shaped dust cloud silhouetted against the bright stellar and nebular field just north of R136 in 
optical images. Again, this cloud is bright in the Spitzer image. It is likely located in the immediate 
foreground of R136. The source IRSW-133 immediately to the south is brighter at J, in which S5 
is invisible (Figure 4a), but the latter is brighter at K and likely dominates in the Spitzer images. 
These sources were also measured at H and Kg by Campbell et al. (2010, their Section 5.4), with 
results similar to ours, although the three available K magnitudes are in relatively poor agreement, 
perhaps due to source confusion or variability. C. Evans (priv. comm.) has pointed out that 
IRSW-133, which Campbell et al. suggested to be a massive star, is VFTS 476 in the current VLT- 
FLAMES Tarantula Survey (Evans et al. 2011); indeed, it has been classified as 0((n)) from the 
presence of He II absorption lines in the VFTS data (Walborn et al. in prep.), with heavy nebular 
emission-line contamination preventing a more precise description. Scale et al. (2009) report Spitzer 
spectroscopy of this source consistent with a YSO nature. 
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3.6. S6 

S6 is the very luminous P2 of Hyland et al. (1992) and NIClGa of Brandner et al. (2001). It is 
a point source in HST optical and IR images, as shown by Walborn et al. (2002) who measured an 
/ mag of 20.9 and placed a V limit of 23. Beware of potential confusion with a Long Period Variable 
(LPV) about 2" NW, which oscillates around 17th mag in / except during dust-formation episodes 
when it becomes fainter than S6 (OGLE-LMC-LPV- 76683; Soszynski et al. 2009). The LPV is 
always much fainter than S6 at longer wavelengths (Figure 4b). To within the HST resolution, S6 
appears to be an instance of monolithic, "isolated" massive star formation. 

3.7. S7 

This intense "double source" is Rubio et al. (1998) IRSN-122/126, the latter dominating Hyland 
et al. (1992) PI; it is also Brandner et al. (2001) NIC12b/d. As extensively discussed by Walborn 
et al. (1999, 2002), they are embedded within a massive dust pillar oriented toward R136 and 
parallel to the adjacent Knot 1 pillar from which a young stellar system has just emerged. Like S4 
and S6, these luminous sources can be measured at / (20.4 for S7A and 21.2 for S7B), with V hmits 
of 22 for both; hence they are also outstanding candidates for future spectroscopic investigation 
over a wide wavelength range. 

3.8. S8 

In striking contrast to S6, this bright source (P4 of Hyland et al. 1992; IRSN-137 of Rubio et al. 
1998; NIClSb of Brandner et al. 2001, but mislabeled there as IRSN-134) is surrounded by a rich 
cluster of much fainter ones as shown by Walborn et al. (2002), indicating qualitatively different 
modes of massive star formation for these two sources. This circumstance creates difficulties with 
the photometry of S8, because the VISTA and IRSF images from which the YJHK magnitudes 
listed in Table la are derived resolve it from the cluster, whereas the Spitzer images do not. Rubio 
et al. distinguish a relatively bright "tail" of the blended cluster toward the SW as the separate 
source IRSN-134, but only the NICMOS images reveal the actual structure of the cluster. We have 
noted that the cluster is surrounded by a thick annulus of emission in the Spitzer /MIPS 24 fim 
image, which warrants further study. 

3.9. S9 

This remarkable source is a constituent of the intricate YSO/ZAMS phenomena associated 
with the Knot 2 complex of Walborn & Blades (1997). With reference to Figure 3 and Section 3.2 
of Walborn et al. (1999): S9 is the brighter of two extended sources aligned across and approxi- 
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mately equidistant from the central, young optical stellar system (Walborn et al. 2002). The major 
axis of S9 has an extent of V\ or 0.25 pc ^ 50,000 AU. Both extended sources display (opposite) 
color gradients perpendicular to their radius vectors from the central system, suggesting (coun- 
terclockwise) rotation. No evidence of the presumed jets impacting the surrounding dust clouds 
to produce them has been found to date. The spectral type of the brightest stellar component in 
Knot 2, VFTS 621 or Parker (1993) 1429, has been found to be 02 V((f*))z from the VFTS data 
(Walborn et al. in prep.), thus a very massive young object. 

It is surprising that S9 would be one of the brightest Spitzer sources in 30 Dor, but then, its 
nature is entirely unknown. Scale et al. (2009) report Spitzer spectroscopy of this source consistent 
with a YSO, although the slit width is at least ^4'' and S9 is not dominant at 8 fim (Fig. 4c). 
Again, the NICMOS images show that it is nonstellar, so that a very large circumstellar envelope 
would be implied in that case, but the point-symmetric counterpart on the opposite side of the 
central stellar system suggests a more exotic interpretation. Moreover, we have found that S9 and 
Knot 2 are apparently associated with the brightest 30 Dor source in all of the Herschel 100-500 //m 
bands, from the HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE; Meixner 
et al. 2010). Finally, there are also two water masers associated with this complex (Ellingsen et al. 
2010), one of which is very near S9, while the other is roughly opposite across the central stellar 
system, approximately mimicking the alignment of the IR "mystery spots" at a lower PA (Figure 4c; 
Walborn et al. 1999). Clearly further detailed, combined analysis of the multiple rich datasets now 
available for this region will prove rewarding. 



3.10. SIO 

At first glance SIO appears to be one of the brightest sources in the Spitzer images, but on 
closer inspection it is somewhat extended, and indeed, the higher resolution VISTA images reveal 
it as three much fainter sources that we designate A, B, C. Moreover, they may comprise part of a 
larger structure associated with the luminous WN star VFTS 682, newly discovered by that survey 
(Bestenlehner et al. 2011). The WN star appears to have excavated a cavity, around the periphery 
of which there is an array of IR sources, including the SIO system. Another bright K source that 
we designate SIOK forms a triangular asterism with VFTS 682, and SIO at the apex (Figure 4c). 
Rubio et al. (1998) IRSN-169 and its fainter NE companion IRSN-170 lie near the opposite side of 
the cavity (midway between VFTS 682 and S9 in Figures 1 and 2); the blended Spitzer data for 
the two sources are ascribed to "Sll" in Table 1 and Fig. 3. There is also a water maser exactly 
coincident with IRSN-169 (Ellingsen et al. 2010). Thus, along with S3/R135, this structure may 
represent another case of a massive WN wind triggering further star formation. 
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4. Discussion 

4.1. Model Fitting and Masses 

We fitted YSO models to the SEDs of selected, unresolved sources. We used the Robitaille 
et al. (2006) grid of 20,000 2D YSO radiative transfer models for masses from 0.1 through 50 Mgun- 
For each model in the grid, SEDs are calculated at 10 viewing angles, resulting in 200,000 pre- 
computed SEDs that can be compared to the observations using the Robitaille et al. (2007) fitting 
tool. Since these models are computed for single objects, we selected those sources from our sample 
that appear single in the near-IR and mid-IR images: S2, S3, S3K, S4, S6, S7A, S7B, and SIOK. For 
the fitting, we used VISTA J and Kg, IRSF i7, and IRAC 3.6-8.0 /xm fluxes as available (Table la), 
as wen as HST /WFPC2 F814W (-/) and F555W (-y) fluxes or upper limits for S4, S6, S7A, 
and S7B (these magnitudes are cited in the individual source discussions above). The fits for these 
sources are shown in Figure 6. All of the YSO fits indicate Class I objects. 

The formal fit for S3K as a star is superior to that as a YSO, so both fits are shown for this 
source. This result is related to the fact that the Spitzer fluxes are smaller than those at shorter 
wavelengths, unlike all the other cases. The stellar fit yields Tgff = 4750 K and Ay = 5.85. While 
the apparent close association of S3K with S3, R135, and nebular structure shown in the figures 
may be a chance alignment, there is also the interesting possibility that it could be a post- YSO, 
PMS object. 

The physical parameters estimated from the fitting are listed in Table 2: best and average 
stellar luminosity, effective temperature, and mass; envelope and disk masses; and interstellar 
extinction. Since the Robitaille et al. (2006) models assume a Galactic gas-to-dust ratio, we scaled 
the results for the envelope and disk masses to the 30 Dor ratio by multiplying them by a factor 
of 2.1 (Dobashi et al. 2008). The averages are over all YSO model fits with P^r data point 
between X^est Xbest + 2' ^^^^ listed in the table. The uncertainties of the averages quoted 
(in parentheses) in Table 2 represent standard deviations of the means. For sources with very high 
X^/pt and only one fitted model (S6 and S7B), there are of course no averages and the derived 
parameters are highly uncertain. For other fits with x^/P^ > 5, the physical parameters may not 
reflect well the actual properties of the sources. Source S3K is also listed in Table 2 as a YSO and 
only its luminosity is provided; however, as already noted this source is more likely to be a star 
with the parameters given above. The Xmin/P^ stellar photosphere fit is ~6, far lower than 

that for the YSO fit, ~109 as listed in the table. 

The models assume single sources; at the distance of the LMC, even the sources that appear 
unresolved in our data could in fact be multiple or even protoclusters (see also Carlson et al. 
2011). Our derived parameters are valid only on the assumption that a single source dominates 
the luminosity at all wavelengths. 
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Fig. 6. — SEDs of eight apparently single IR sources. The solid black line shows the best-fitting 
model. The dashed line corresponds to the (extincted) SED of the stellar photosphere in the 
best-fitting model. S3K has been fitted as both a YSO and a star, with the latter fit appearing 
preferable. The plots are labeled with Robitaille et al. (2006) model numbers and global values. 
Further parameters resulting from the fits are given in Table 2. 



Table 2. The SED Fitting Results: Physical Parameters 



Source U (10^ Lq) (10^ K) (Mq) Menv (10^ Mq) Maisk (Mq) Ay X^in/Pt^ Hfits" 

Name best ave best ave best ave best ave best ave best ave 



52 3.4 7.7(0.7)^ 1.2 2.1(0.1) 20.1 23.0(0.7) 2.2 6.0(0.7) 0.15(0.05) 6.5 5.4(2.7) 0.4 109 

53 8.2 10.5(2.5) 3.8 3.4(0.6) 25.2 26.7(3.2) 13.3 14.4(7.7) 10.0 9.7(0.4) 14.7 3 

S3K .... 1.7 ••• ••• ••• ••• ••• ••• ••• ••• ••• ••• ••• 109.3 1 

54 10.7 8.2(2.5) 3.9 3.8(0.1) 27.4 24.5(2.9) 10.9 6.4(4.5) 0.19 0.10(0.09) 1.8 1.2(0.6) 8.7 2 

S6 3.7 3.4 ••• 18.4 ••• 12.0 ••• 0.13 ••• 3.0 ••• 37.9 1 

S7A .... 3.0 2.7(0.4) 2.0 2.4(0.4) 15.2 16.4(0.6) 1.6 3.1(1.2) 0.04 0.13(0.09) 3.2 2.8(1.7) 1.4 9 

S7B .... 7.3 ••• 1.7 ••• 23.4 ••• 71.5 ••• ••• 0.1 ••• 41.1 1 

SIOK ... 0.7 2.4(1.8) 2.5 2.0(0.1) 9.1 10.9(1.5) 0.5 2.8(1.9) 0.28 0.11(0.02) 4.4 4.5(2.4) 0.5 24 



^ The uncertainties represent standard deviations of the means. 

^ Xmin /pt is a P^r data point for the best-fit model. 

^ Ufits is the number of fits with X'^/P^ between Xmin /P^ Xmin /pt+2. 
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5 . S ummar y / O ut lo ok 

Based primarily upon Spitzer /SAG^, VISTA/VMC, and IRSF data, with reference to HST / 
NICMOS/WFPC2/WFC3 images for higher spatial resolution, we have presented the first compre- 
hensive description of the brightest IR point sources throughout the triggered star formation within 
the 30 Doradus Nebula in the LMC. The Spitzer data extend prior work to longer wavelengths, and 
all of the subsequent IR images cover fields not included in the NICMOS dataset, thus completing 
the global IR view of 30 Dor. The overall spatial distribution of these sources is consistent with 
the two-stage starburst picture discussed previously, but the intricacy and diversity among the in- 
dividual "point" sources revealed here is striking. There are clear implications for all YSO studies 
in the Magellanic Clouds, not to mention at greater distances. In summary, the most outstanding 
individual source results are listed below: 

• SI is a newly resolved compact H II region containing a multiple system of massive stars 
readily seen in the / band, possibly triggered by adjacent SN events. The object is named 
"The Skull Nebula" from its visual appearance in the WFC3 image. 

• S3, the multiple SIO, and other nearby sources are apparently associated with isolated, lumi- 
nous WN stars and may be the result of local triggering by the massive winds of the latter. 
A possible analogous instance in the Galaxy was recently reported by Liu et al. (2012). 

• S4, S7A, and STB are located within massive dust pillars oriented toward R136, as discussed 
in previous work. 

• S6 is a very luminous, unresolved point source in the HST images and may represent an 
instance of isolated massive star formation. 

• S8 is the most luminous source in a small cluster of fainter objects, in marked contrast to S6. 
The cluster is surrounded by a 24 fim shell in the Spitzer /MIPS image. 

• S9 is an unprecedented extended source of unknown nature; with its fainter counterpart on the 
opposite side of a very young, massive optical system, its morphological appearance suggests 
the impact point of an invisible, rotating jet on surrounding dust clouds. Surprisingly, it 
is one of the brightest Spitzer sources, and it is related to the brightest Herschel source in 
the region. There are also two water masers among the array of interesting objects in this 
configuration, which will likely provide unique information about massive star formation. 

We have fitted YSO models to the data for eight apparently (predominantly) single objects to 
derive physical parameters including stellar masses, which fall in the 10-30 Mq range, in general 
agreement with the results of Hyland et al. (1992) for four of the brightest sources. 

These results can be extended in several ways, at least two of which are obvious. First, 
measurements and analysis of the complete IR sample in this region should be carried out in order 
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to estimate the total mass of the second, triggered generation in 30 Doradus, for comparison with 
the first generation to illuminate the general structure and evolution of starbursts on this scale. 

Second, broad- wavelength spectroscopy of the bright sources is essential to advance understand- 
ing of their diverse natures. A current ESO/VLT/X-Shooter program (PI L. Kaper) is addressing 
this objective. As noted in the individual discussions, some of these sources are detectable at / and 
even at which promises significant further insights into their physical properties. For example, 
a comparable observation of a Galactic YSO has recently yielded photospheric absorption as well 
as circumstellar emission features from the optical through the NIR (Kaper et al. 2011; Ochsendorf 
et al. 2011). 

This paper is based largely on (SAGE) observations made with the Spitzer Space Telescope^ 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under 
contract with NASA; and in part on ( VISTA/ VMC) data products from observations made with 
ESO Telescopes at the La Silla or Paranal Observatories under ESO program ID 179.B-2003; 
and in part on data from the IRSF Magellanic Clouds Point Source Catalogue obtained at the 
SAAO. Ancillary use has been made of data products from 2MASS, which is a joint project of the 
University of Massachusetts and the Infrared Processing and Analysis Center/California Institute 
of Technology, funded by NASA and NSF. Thanks to Elena Sabbi and Danny Lennon for their 
unpublished HST /WFC3 data and discussion regarding SI. M.S. acknowledges financial support 
from the NASA ADAP award NNX11AG50G. Publication support was provided by NASA through 
grant GO-12465.01 (PI P. Crowther) from STScI. 
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